NBI-HE-95-05 
February 1995 



Scaling in four dimensional quantum gravity 

J. Ambj0rn and J. Jurkiewicz 1 

The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



Abstract 



We discuss scaling relations in four dimensional simplicial quantum gravity. Using nu- 
merical results obtained with a new algorithm called "baby universe surgery" we study 
the critical region of the theory. The position of the phase transition is given with high 
accuracy and some critical exponents are measured. Their values prove that the transi- 
tion is continuous. We discuss the properties of two distinct phases of the theory. For 
large values of the bare gravitational coupling constant the internal Hausdorff dimension 
is two (the elongated phase), and the continuum theory is that of so called branched poly- 
mers. For small values of the bare gravitational coupling constant the internal Hausdorff 
dimension seems to be infinite (the crumpled phase). We conjecture that this phase corre- 
sponds to a theory of topological gravity. At the transition point the Hausdorff dimension 
might be finite and larger than two. This transition point is a potential candidate for a 
non-perturbative theory of quantum gravity. 
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1 Introduction 



Many aspects of two-dimensional quantum gravity are by now well understood. Using 
Liouville theory [1] we can understand how to quantize two-dimensional quantum gravity 
coupled to conformal field theories with a central charge c < 1. In addition it is understood 
how to discretize the theory in a way consistent with the underlying reparametrization 
invariance of the continuum theory [2-5]. The continuum theory is recovered at the critical 
point of discretized theory. In this way it was possible to utilize the connection between 
critical phenomena in solid state physics and quantum field theory, which previously has 
proven so useful. The unique situation if c < 1 is that a number of the discretized 
theories allow an explicit solution. In this way it is possible to compare in detail the exact 
results with the results of numerical simulations of the discretized systems. Monte Carlo 
simulations have always been an important tool in the theory of critical phenomena and 
it can be verified that such methods work very well for two-dimensional quantum gravity 
coupled to conformal matter too (for a recent review, see [6]). This observation allows 
us to use numerical simulations in situations where two-dimensional gravity coupled to 
matter cannot be solved analytically. The numerical approach provided us with a good 
understanding of the conformal field theories with c ^> 1 coupled to quantum gravity [7], 
and it served as an important inspiration for analytic approximations in this region [8-11]. 

The discretized version of two-dimensional quantum gravity mentioned above is known 
as the theory of dynamical triangulations or simplicial quantum gravity. It can be general- 
ized to higher dimensions (in the three-dimensional case see [12,13], in the four-dimensional 
case [14-24]). By such a generalization we enter unchartered territory, since we have no 
continuum theory of Euclidean quantum gravity with which we can compare. Nevertheless 
the generalization is so simple and aesthetically appealing that it deserves a serious inves- 
tigation. It has not yet been possible to use the strong analytic methods available in two 
dimensions but as we show in this article, many aspects of the theory can be investigated 
to a sufficient degree of precision by computer simulations. 

The rest or the article is organized as follows: In section 2 we review the definition of 
the model. In section 3 we discuss the general scaling relations which should be fulfilled if 
a continuum limit exists in the conventional setting of critical phenomena. These scaling 
relations are necessary in order to understand the critical behavior of the model. In 
section 4 we report on the extensive numerical results obtained so far. Subsection 4.1 is a 
description of a new algorithm for Monte Carlo simulation. It is very efficient in so-called 
elongated phase where it eliminates the critical slowing down and allows us to determine 
critical exponents with great precision. Finally section 5 contains a discussion of results 
and the future perspectives. 

2 The model 

The idea behind simplicial quantum gravity in more than two dimensions goes back to 
Weingarten [25] and found its "modern" form in [12] in three dimensions and in [14, 15] in 
four dimensions 2 . In the continuum formulation the task is to perform the path integral 

2 See [16, 18] for reviews. 
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over equivalence classes of metrics with the action given by that of general relativity: 



Z(A,G,X,h)= V / V ^ P -5 g [A,G,S] 



^ )m{y.) Vol(diff) 



x(s)eTo P 

5a [G, A, S] = / rf 4 ?^ ( A ) + boundary terms. (2) 

Jm(E) V 167rG / 

Here S denotes a three-dimensional boundary and the boundary terms depend on the 
induced metric h^ v on the boundary and the extrinsic curvature induced on S. It is 
needed in order that the action satisfies the equation: 

+ Sg 2 = Sg 1+ g 2 (3) 

which is a shorthand notation for the following: The sum of the action for a manifold with 
metric gi and boundary S and induced metric h and the action of a different manifold 
with metric g 2 and opposite oriented boundary — S and induced metric h must be equal 
to the action of the manifold constructed by gluing together the two manifolds along S. 
The condition (3) is needed in order that the partition function Z(T<i, hi; T, 2 , h 2 ) satisfies 
the fundamental composition law: 

hi, S 3 , h 3 ) = ^2 Z(T, 1 ,h 1 ;T, 2 ,h 2 ) Z(T, 2 ,h2;T,3,h 3 ), (4) 

which expresses the fact that we should be able to sum over intermediate states, charac- 
terized here by the boundary T, 2 and the equivalence class of three-metrics h 2 (T, 2 ). 

In (2) we are instructed to perform the summation over all manifolds with fixed bound- 
aries and for each manifold to perform the integration over equivalence classes of Rieman- 
nian metrics. Nobody has yet been able to make sense of this expression. There is a 
number of reasons for this. Firstly it is not clear what is meant by the summation J^Top 
over topologies. In four dimensions there is no equivalence between smooth structures and 
topological structures. There even exist compact topological manifolds with a countable 
infinity of non-equivalent smooth structures and there exist topological manifolds which 
do not admit any smooth structure. If we restrict ourselves to smooth structures we have 
the additional problem that the Euclidean action is unbounded from below. This implies 
that it is not even clear how to define the path integral over equivalence classes of metrics 
in the case of a fixed structure. Even if we could ignore the problem of the unboundedness 
of the Euclidean Einstein action we would still be faced with the problem of evaluation 
of the path integral. Apart from the technical difficulties of dealing with reparametriza- 
tion invariant theories we also face the unpleasant fact that the Einstein action is not 
renormalizable. 

Clearly these problems await the proper non-perturbative formulation, either by being 
embedded in a larger theory (like string theory), or by some discretization method which 
can tame (1). Simplicial gravity is a suggestion of a discretization which allows us to 
discuss some aspects of (1) in an non-perturbative setting. It attempts to approximate 
the integration over equivalence classes by the summation over combinatorially equivalent 
piecewise linear manifolds. 
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Here the following should be noted: in four dimensions the equivalence classes of 
piecewise linear structures are in one-to-one correspondence with smooth structures. If 
by J2top we understand the summation over all smooth structures, we change nothing by 
replacing "smooth" with "piecewise linear". To each abstract piecewise linear structure 
we can associate a metric in the following way: Consider each link in the triangulation to 
be of equal length a (the "lattice spacing"). If we consider the interior of each four simplex 
as flat we have defined the concept of distance on the piecewise linear manifold, and in 
addition Regge calculus associates in a natural way an action to such a piecewise linear 
manifold. In this case the expression for the action becomes very simple. The volume is 
proportional to the number of four-simplexes and the curvature term, which is the sum of 
deficit angles of the two-simplexes, will be proportional to the number of four-simplexes 
containing the two-simplex. For a closed manifold the sum over all deficit angles can 
now be expressed as a linear combination of the total number of four-simplexes, N 4 , and 
the total number of two-simplexes, N 2 . If we take the lattice unit a = 1 the continuum 
action (2) for a piecewise linear manifold characterized by an abstract triangulation T 
with boundary B takes the exceedingly simple form: 

St[K h] = k 4 N 4 {T) - k 2 N 2 {T) + h 2 N 2 {B), (5) 

where k 2 ~ 1/G. One could ask if such a simple expression has any chance of representing 
a theory of quantum gravity. Our attitude is the opposite: It is intriguing that there exists 
such a simple expression and hopefully the simplicity might lead to the explicit solution 
of quantum gravity by analytic means in the future. This is precisely what has happened 
in two dimensional quantum gravity. 
The path integral (1) is replaced by 

Z[k 2 ,k 4 ] = J2 7 ^e- ST[k2 > k4 \ (6) 

where Cy is a symmetry factor of the graph T (the order of the automorphism group of 
the graph). The integration over equivalence classes of metrics in this formula has been 
replaced by the summation over a grid in the space of equivalence classes of metrics. Each 
abstract triangulation (up to graph isomorphism) represents such a point, the metric 
being assigned as described above. Different triangulations (up to graph isomorphism) 
correspond to inequivalent metrics and consequently to different points. The weight of the 
given metric (or triangulation) is given by the exponential of the action (the Boltzmann 
weight), calculated by Regge's prescription. The basic assumption is that this grid (with 
Boltzmann weight included) becomes dense in the total space of equivalence classes of 
metrics (with Boltzmann weight included) for certain values of the bare coupling constants 
k 2 ,k 4 , such that the sum (6) provides us with a realization of (1). 

One annoying aspect of (1) is the split in a summation over topologies and the inte- 
gration over Riemannian structures. It seems that in a natural way we have been able 
to combine these in simplicial quantum gravity by formula (6). This is yet another fea- 
ture of simplicial quantum gravity which has been very important in the two-dimensional 
models, where the summation over topologies, known as the double scaling limit, can be 
performed in the simplicial quantum gravity models. Unfortunately the sum (6) is still 
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quite formal, and in more that two dimensions it is not yet known how to make sense of it. 
If we insist on an unrestricted sum over topologies there is a huge entropy factor, which is 
not sufficiently damped by the action, no matter what are the values of the bare coupling 
constants k 2 , k 4 . In this paper we shall make no attempt to define the sum over topologies, 
but confine ourself to the simplest of all cases, that of four-manifolds of spherical topology, 
except for simple boundaries each of a topology corresponding to S* 3 . 

At the moment we restrict the topology in the summation (6) the sum becomes well 
defined in some range of the (k2,k 4 ) coupling constant plane. This follows from the fact 
that the number of triangulations of a given topology seems to be exponentially bounded 
with the volume, i.e. the number of four-simplexes. The analytic proof is well known in 
two dimensions [26]. In higher dimensions recently there has been some controversy about 
this point [27]. Computer simulations support the exponential bound [28,21] and very 
recently there appeared an analytic proof of the existence of the exponential bound in 
all dimensions and for all (fixed) topologies [31]. The existence of the exponential bound 
implies that the number of triangulations A/" (^4) for the number of four-simplexes equal 
to N 4 (and for fixed topology) will be bounded by 

N(N 4 ) < e cNi . (7) 

Since we have a trivial bound N 2 (T) < constant- N 4 (T) it follows that for any fixed k 2 (6) 
is well defined for sufficiently large k 4 . For a fixed k 2 there is a critical lower value k1(k 2 ) 
of k 4 such that the partition function (6) is well defined for k 4 > k 4 (k 2 ) and divergent for 
k 4 < k1(k 2 ). To be more specific we can introduce the partition function for fixed (lattice) 
volume: 

Z(k 2l N 4 ) = £ e k * N ^ = f k2 (N 4 )e k ^ N \ (8) 
TeT(JV 4 ) 

where T(N 4 ) denotes the class of triangulations of fixed topology and boundary conditions, 
here always taken to be that of a four-sphere with S* 3 boundaries. The function fk 2 (N 4 ) 
denotes a subleading factor and it will be discussed in detail later. Here it is sufficient to 
say that it will depend on the boundary conditions and on k 2 . 

In the limit k 4 — > k 4 (k 2 ) large N 4 will dominate (6). It corresponds to the infinite 
(lattice) volume limit. This implies that the infinite volume limit is found by the approach 
to a the critical line k 4 = k 4 (k 2 ) m the (k 2 ,k 4 ) coupling constant plane. It should be 
emphasized that the existence of an infinite volume limit does not necessarily imply that 
we have an interesting continuum limit. Rather, the situation should be compared with a 
lattice theory. For a spin system we can always take the infinite volume limit by taking 
the lattice volume to infinity. However, only for a specific range of the coupling constants 
we shall be able to obtain a spin-spin correlation length which diverges with respect to 
the lattice spacing, and it is only for this range of coupling constants that we have any 
hope of obtaining a conventional continuum field theory. 

3 Scaling relations 

In the case of a spin system on a regular lattice the spin-spin correlation length, or the 
mass gap m(P), is the quantity which characterizes the scaling behavior. If m(/3) — > for 
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fi —7- f3 c , where fi c denotes the critical point, it might be possible to define a continuum 
limit of the spin system, since the correlation length becomes infinite compared to the 
lattice spacing. All scaling relations can be derived from this scaling hypothesis. We can 
introduce a similar object in simplicial quantum gravity. 

Let us first discuss the formal continuum physics of the object. In order to address the 
question of divergent correlation length we have first to introduce the concept of length 
in quantum gravity. It can be done even if we have to integrate over the metrics. The 
simplest example of a reparametrization invariant correlation function which depends on 
a geodesic distance is the volume-volume correlator: 

G{r; A, G) = (// \fg^)d^ yjg&d*? 8(d g (Z,?) - r)} , (9) 



wf 



here d g (£, £) denote the geodesic distance between the points £ and calculated with the 



metric g^ v , and where the average is performed with the formal path integral of quantum 
gravity, i.e. over all equivalence classes of metrics [g]. The subscript c in (-) c denotes 
the connected part of the correlator. In a similar way one can include in the average 
reparametrization invariant local operators like the scalar curvature R(x) or B^ VK p - 

Gn{r; A, G) = (// y/gi^R(0^ yfi^R^d 4 ^ 8(d g (Z,?) - r)} . (10) 

The integrals of G(r, A,G) and G#(r; A, G) are directly related to the partition function 
Z(A, G) itself by: 

°° a rt a r\ d 2 log^(A,G) 

dr G(r;A,G) = — 2 (11) 

and 

In principle one could also consider such correlators for a fixed volume V of the universe, 
e.g. 

G(r;V,G)=(J f yfii^d 4 ^ yfii^jd 4 ? 8(d g (Z, £') - r)5(f d 4 ^ - V)} . (13) 

While not so interesting from a physical point of view the knowledge of the correlators for 
any volume allows us to reconstruct the former correlators by e.g.: 

G(r;A,G)= / dVe~ AV G{r; V, G), G{r;V,G)= / - — : e AV G{r; A, G). (14) 

JO Jc-ioo 27T2 

From a numerical point of view G(r;V,G) are more convenient objects. 

The above mentioned correlation functions have a direct and simple translation to the 
discretized approximation where all quantities are well defined. Consider piecewise linear 
manifolds with two marked four-simplexes, or equivalently with two boundaries which are 
the boundaries of two four-simplexes. The geodesic distance between any two points is 
uniquely defined from the piecewise linear structure and the requirement that the link 
length is a = 1. Let us for simplicity (in particular for the numerical work to be reported 
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Figure 1: Graphical illustration of inequality (16) in the two dimensional case. Eq. (16) 
corresponds to the situation where the boundaries are contracted to points, but is in fact 
valid for arbitrary boundaries. 



below) adopt a simplified notation 3 : We define a path between two four-simplexes as a 
sequence of neighboring four-simplexes connecting the two. The geodesic distance is the 
shortest such path. This definition can in an obvious way be used to define the geodesic 
distance between two boundaries. Let us introduce the notation ds^ for the boundary of 
a four-simplex S4. We can now define the 2-point function as the partition function (1): 

G(r; k 2 , k 4 ) = ^ e-Mfc,*^ (15) 

TGT 2 (r) 

where T 2 (r) denotes the class of triangulations of S* 4 with two boundaries ds4(l) and 
ds4(2) separated by a geodesic distance r. The important property of G(r; k 2 , ^4) is its 
exponential decay with the geodesic distance r. It follows from the inequality: 

G(ri + r 2 ; k 2 , k 4 ) > const • G(ri, k 2 , k 4 ) G(r 2 ; k 2 , k 4 ). (16) 

This inequality has a simple geometrical interpretation. Given two triangulations T\ with 
boundaries ds±(l) and ds^^i) separated by distance r\, and T 2 with boundaries ds4(4) 
and ds4(2) separated by distance r 2 , we can clearly construct a triangulation LL T\ + T 2 
with boundaries ds±(\) and 9s4(2) separated by a distance r\ + r 2 by gluing together T\ 
and T 2 along 9s4(3) and d 84(A). Due to the property (3) the action for the LL T\ + T 2 V will 
be the sum of the actions for T\ and T 2 . Except for a trivial constant factor of order one 
due to different ways of gluing ds4(3) and ds4(4) together each T\ and T 2 entering on the 
lhs of eq. (17) is present in the sum on the rhs of (16) in the form LL T\ + T 2 V with the same 
Boltzmann weight. In addition there are of course many terms in the sum on the rhs of 
(16) which have no representation on the lhs of (16). The inequality is shown graphically 
in fig. 1. 

The important implication of (16) is that the function — logG(r; k 2 , ^4) is a subadditive 
function of r. This implies that 

lim — — — — — — = m(k 2 , k^) exists. (17) 

r— >oo f 

3 This definition of geodesic distance has the additionai advantage that it appiies to more generai com- 
plexes than genuine triangulations. In the broader context of matrix models one will be forced to sum 
over such general complexes [12] and here it is convenient to be able to use the same concept of geodesic 
distance. 



7 



In addition it is not difficult to show by combinatorial arguments that for k 4 > k 4 (k 2 ) 

> c 2 G(r;k 2 ,k 4 ), (18) 



dG(r = l;k 2 ,k 4 ) 2 



dk 4 

dG(r;k 2 ,k 4 ) „ 2 
Wk 4 G(r;k 2 ,k 4 ) (19) 

where c 2 and c 2 are two constants of O(l). The arguments are identical to the ones 
originally used for random surfaces [32] and they imply that 

m(k 2 ,k 4 )>0 and gm ( fc 2, k 4 ) > Q ^ k 4 > k c 4 {k 2 ). (20) 

O K 4 

The first inequality follows from (17) and (18) in the limit r — > 00. The second inequality 
for from (17) and (19) for r — > 00. 

To summarize: For k 4 > k1(k 2 ) the 2-point function falls off exponentially, but although 
the mass m(k 2 , k 4 ) is decreasing for k 4 — > k 4 (k 2 ) the general arguments do not tell us that 
the mass scales to zero. 

Let us in the following assume that 

m(k 2 , k 4 ) ~ {k 4 - k c 4 {k 2 )Y for k 4 k c 4 {k 2 ), (21) 

and deduce the consequences of this assumption. In the following we will suppress the 
explicit k 2 dependence and introduce the shorthand notation Ak 4 = k 4 — k 4 (k 2 ). Let us 
further introduce the 2-point function G(r; N 4 ) for a fixed volume N 4 : 

G{r;N 4 )= z k2N2{T) - (22) 

TeT 2 (r;N 4 ) 

The function is convenient in numerical simulations as will become clear in the following. 
It is related to G(r; k 4 ) by a (discrete) Laplace transform: 

G(r-k 4 ) = ^- kiNi G(r-N 4 ). (23) 

From the large distance behavior of G(r, k 4 ) dictated by (17) and (21) we deduce the long 
distance behavior of G(r;N 4 ) by an inverse Laplace transform of (22): 

G{ n k 4 )~e- cr ^ v => G{ n N 4 )^e-'^l N ^ l{1 - v) (24) 

The mass exponent v has the interpretation as the inverse internal Hausdorff dimension 
of the ensemble of four-manifolds: 

v=l/d H . (25) 

Let us prove this relation. Given the ensemble of four-manifolds with boundaries ds 4 (l) 
and ds 4 {2) separated by a geodesic distance r it is natural to define the Hausdorff dimen- 
sion by 

(N 4 ) r ~ r dH (26) 
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here (A(N4)) r is defined by 



(A(N )) ~ ^eT 2(r) A(N 4 (T))e-^ 

{A{N4))r ~ ET £ r 2(r) e-^W (2?) 

For large r we can use the asymptotic form G(r; k 4 ) ~ e _m ( fc4 ) r of the 2-point function to 
calculate the (N 4 ): 

This formula expresses the trivial fact that for fixed > k\ and r — > 00 the typical 
manifold is a long tube and the Hausdorff dimension is one. The same would be the case 
for an ordinary random walk and it highlights the way the scaling limit should be taken: 
We want r — > 00 in order to be able to ignore lattice artifacts, but at the same time we 
have to scale m(k 4 ) — > such that the 2-point function ~ e _m ( fc4 ) r survives, i.e.: 

m(k4)r = const. or (Ak 4 ) u r = const. (29) 

If we use this in (28) we get: 

<Ar 4 > r ~ r 1 /", i.e. d H = k (30) 

One can compare this result with the ordinary random walk in R n where it is well known 
that the mass exponent v =1/2 and the {external) Hausdorff dimension is two. 

While G(r; k 4 ) is a convenient object from a theoretical point of view, it is G(r; N4) 
which is readily available from the point of view of computer simulations. Let us assume 
the Hausdorff dimension of our ensemble of manifolds is dfj- An alternative definition of 
dfj would be the following: Let us consider the ensemble of four-manifolds with a fixed 
volume N4 and one boundary ds^. Let us move away from the boundary in "spherical" 
shells of thickness 1 and count the number of four-simplexes n(r) in the shells a distance 
r from the boundary [14,23]. If we take the average over the ensemble of manifolds we 
expect the following behavior, as long as r is not too large compared to N^ dH : 

(n(r)) N4 ~ r^" 1 for 1< r < N l J dH . (31) 

The number n(r) is closely related to G(r;N4). Since we can view G(r = 0; N4) as the 
partition function for the ensemble of manifolds with one boundary 834 except for a trivial 
combinatorial factor of order one, we have: 



G(r;N 4 
G(0; N 4 ) ' 



Wr)) Ni ~ ^f^, (32) 



from which by (24) we deduce the long distance behavior of {n(r)) N : 

-(rN lld "Y Hl{dH ~ 1) ,„ 

(n(r)) N ~ e V 4 J for N l J dH < r < N 4 . (33) 
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While the long distance behavior of G(r; N 4 ) is determined from that of G(r;£; 4 ), which 
again follows from the scaling assumption (21), the short distance behavior of G(r; N 4 ) is 
dictated by (31): 

G(r; N 4 ) ~ G(0; N 4 ) r^" 1 for 1< r < N l J dH , (34) 

and we can determine the short distance behavior of G(r; k 4 ) from (23). Recall the ex- 
pression (8) for the partition function Z(N 4 ) for finite volume (we still suppress explicit 
reference to A^). For N 4 n °t too small we have 

G(0;7V 4 ) ~ N 4 Z{N 4 ), (35) 

since G(0, N 4 ) corresponds to manifolds with one boundary ds 4 and this boundary can be 
inserted at any of the N 4 simplexes. For very small symmetric triangulations additional 
combinatorial factors might be present, this is why we say that N 4 should not be too 
small. According to (8) we expect Z(N 4 ) to grow exponentially with N4. The critical 
behavior of G(r, k 4 ) is related to the subleading behavior f(N 4 ). Two functional forms of 
the subleading behavior have been suggested: 

f(N 4 ) = (36) 

and 

f{N4)=e -cN^ (3?) 

The critical behavior corresponding to these two possibilities will be very different. 

Let us first discuss the implications of the form (36) which is well known from the 
string theory and two-dimensional simplicial gravity. Like in two dimensions [9, 8] one 
can prove that either 7 < or 7 = 1/n, n > 2. For r <C I/A&4 we get, using (23) and 
(34)-(36): 

-(rN lld A dHl{dH ~ l) 

Gfr.y-A-^iV^e-^e i 4 ) -r^" 1 . (38) 

It follows that the short distance behavior of G(r, A; 4 ) is determined from the Hausdorff 
dimension and the "entropy exponent" 4 7. This behavior reflects both kinds of fractal 
structures readily definable in quantum gravity: The first is the internal Hausdorff dimen- 
sion dfj and the second is the baby universe structure. 

The fact that the exponent 7 defined by (36) reflects the distribution of baby universes 
was shown in [33] in two dimensional simplicial gravity, and generalized to four dimensional 
simplicial quantum gravity in [19]. Let us first define a minimal bottleneck baby universe, 
abbreviated "minbu" . A "bottleneck" in a triangulation of S* 4 is a collection of five three- 
simplexes isomorphic to the boundary 9s 4 of a four-simplex. If we cut our piecewise 
linear manifold along such a closed three dimensional manifold we will separate it in two. 
Usually this separation will be trivial in the sense that the three-dimensional manifold 

4 The notation "entropy exponent" is justified in case where &2 = 0, since Z(Ni) in this case simply 
counts the number of triangulations. 
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will be the boundary of a four-simplex S4 in the triangulation and all we have done is 
to separate this four-simplex from the rest of the triangulation. However, in case the 
three-dimensional closed manifold is not the boundary of a four-simplex belonging to the 
triangulation, cutting along the boundary will produce a separation into two non-trivial 
parts of the triangulation. The smaller one is called a "minimal bottleneck baby universe", 
abbreviated "minbu" as mentioned above, the larger the "parent". 

The average number of minbu 's of volume V 4 on 4-manifolds of total volume N4 can be 
expressed in terms of the partition function Z(N) and will be given by a general formula 

120 

(Af(V 4 )) Ni « ^^V 4 Z(V 4 ) (N 4 - V 4 )Z(N 4 - V 4 ) (39) 

where 120 is the number of ways one can glue the two boundaries of the minbu and the 
parent together with opposite orientation of the boundary. The additional factors V 4 and 
N 4 — V 4 appear because the minbu and the parent both have a minimal boundary ds 4 , 
i.e. for a generic large surface of volume V 4 and no accidental symmetry factors there will 
be V 4 such manifolds for each closed manifold since the boundary can be placed at any of 
the V 4 four-simplexes. If we assume that the canonical partition function is given by (36), 
we get: 

(Af(V 4 )) Ni ~ N 4 [V 4 (l - V4/N4)] 1 - 2 . (40) 

which has the obvious interpretation that the probability of having a branching to a minbu 
of volume V 4 per unit volume is V^~ 2 (l — V 4 /N 4 )' y ~ 2 . As shown in [34,19] "minbu counting" 
is a very convenient tool for measuring 7 if 7 is not too negative, both in two-dimensional 
and four-dimensional simplicial quantum gravity. 

The "entropy exponent" 7 is also called the "susceptibility exponent". The reason is 
that it appears in the integral (11) of the two-point function G(r;k4): 

^G(r; k 4 ) = xv(k 4 ) = \c\-y + less sm gul ar terms. (41) 

r V 4 4/ 

From the continuum formula (11) it follows that 

, d 2 log Z(k 4 ) . . 

Xv(k4) = |^ (42) 

It follows that the volume- volume fluctuations diverge if 7 > and stay finite if 7 < 0. 

It is also possible to define the discretized version of the continuum curvature-curvature 
correlation function (10). The integral of this discretized curvature-curvature correlation 
function x(^2) nas been extensively used in the former studies of possible phase transitions 
as a function of k 2 . The analogue to the continuum formula (12) would be: 

Mh) = d2l ° g d * {k2 \ k 4 = kt(k 2 ), (43) 

and the existence of a k\ where XR^kq) diverges would be an indication of a phase transi- 
tion. In this study we will not focus on this quantity, but can still add to the understanding 
of its the behavior (see section 4.6). 
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The asymptotic form (37) leads to a scaling behavior different from the one just discussed. 
Formula (39) is of course valid also in this case. In the following we shall see that there 
is good evidence that the behavior (37) is related to dfj = oo (or v — > 0). From (24) we 
get: 

G(r; N 4 ) -> f(N 4 ) e - Ak * N *- cr for d H -> oo, r > r (44) 

where ro is some short distance cut-off of 0(1) with only a weak N 4 dependence. The r 
dependence factorizes and results in the lack of scaling of the mass: 

G(r;k 4 ) = Y^G(r;N 4 ) e~ k * N * ~ e" cr . (45) 

In principle it is still possible to maintain a scaling of the mass m(k 4 ) for k 4 — > k 4 even if 
the Hausdorff dimension is infinite: A behavior like 

m{k 4 ) ~ 1 x (46) 
lo S Afc7 

is possible. However, it will not lead to a factorization of the r dependence in G(r;N 4 ), 
but to a dependence: 

r i r^l lo S lo s( N il T ) \ 

G(r;7V 4 )~e 1 °s(n 4 /t) i°g(w 4 /r) >_ (47) 

Although this is a weak N 4 dependence we shall be able to distinguish between (45) and 
(46) in the numerical simulations. 

It is clear that it is meaningless to talk about a susceptibility exponent 7 in case we 
have the asymptotic behavior (37). Formally 7 = —00. While a negative exponent 7 
implies that the susceptibility Xv{k 4 ) of volume fluctuations remains finite at the critical 
point k 4 , the derivative Xy(^4) after k 4 will diverge for sufficiently large n. This will not 
be the case if we have the asymptotic behavior (37). 



4 The computer simulations 

Since numerical simulations will be the source of information about the properties of the 
four-dimensional simplicial systems let us start this section with a discussion of certain 
theoretical and practical problems connected with Monte Carlo simulations of four di- 
mensional quantum gravity. Given an initial triangulation To(Ai) of a four-manifold Ai 
the Monte Carlo simulations described here make use of two kinds of changes, a standard 
finite set of so called local moves (there are five different local moves for triangulations 
of four-manifolds) and global "baby universe surgery" moves described below. The lo- 
cal moves were discussed in a number of articles 5 Each of them can be implemented by 
a finite number of computer instructions, independently of the size of the triangulation 
and their weights are chosen as the Boltzmann weights dictated by the classical action of 
general relativity, i.e. by e~ Sr . The local moves are ergodic. By this we mean that given 
any two triangulations of the manifold it is possible in a finite number of moves, carried 
out successively, to change one of the triangulations to the other. The local moves seem 
therefore to provide a perfect setup from a computational point of view. 

5 For reviews see e.g. [19, 18]. 
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There is however one serious theoretical problem related to the fact that there exist 
four-manifolds which are algorithmically unrecognizable. Denote such a manifold by Aio 
and let it be finitely presented by a triangulation T(Aio). The algorithmic unrecogniz- 
ability of .Mo means that there exists no algorithm which allows us to decide whether 
another manifold Ai, again finitely presented by a triangulation T(Ai) is combinatorially 
equivalent 6 to .Mo- When this is combined with the existence of the finite set of local 
moves which is able to bring us between any two triangulations of Aio in a finite number 
of steps, but where this number is a function of the chosen triangulations, one can prove 
the following theorem [35]: 

The number of moves needed to connect two triangulations of Aio, T and T' with N^T) = 
N^T'), cannot be bounded by any recursive function r(N4). 

This theorem implies that there will be very large barriers between some classes of tri- 
angulations of Aio and there would be triangulations which can never be reached in any 
reasonable number of steps. Of course the number of configurations which are separated 
from some standard triangulation of Aio by such barriers could vanish relative to the total 
number of configurations as a function of N4. In ref. [35] it was conjectured that it will 
not be the case. If the conjecture is correct a Monte Carlo method based on the finite 
set of local moves will never get around effectively in the class of triangulations of .Mo- 
We can say that the moves, although ergodic in the class of triangulations, will not be 
computationally ergodic [35]. 

The situation for S* 4 which is the manifold we use in the actual Monte Carlo simula- 
tions was discussed in [28]. It is unknown whether S* 4 is algorithmically recognizable in the 
class of four-manifolds. If S* 4 is algorithmically unrecognizable the arguments given above 
for Aio apply. Any attempt to see high barriers separating triangulations for S* 4 has failed 
so far 7 . Does this mean that we have numerical evidence that S* 4 is algorithmically rec- 
ognizable? Unfortunately not. While the four-dimensional manifolds which are presently 
known to be algorithmically unrecognizable are rather unwieldy, the situation changes if 
we move one dimension up, and even a simple manifold like S 5 is algorithmically unrec- 
ognizable. Computer simulations for S 5 give results similar to those of S* 4 , i.e. up to now 
there has been no trace of any barriers [36]. Two interpretations are possible: Either the 
class of configurations of S 5 which are separated from the rest by high barriers contains 
very special configurations which are of measure zero and not important for numerical 
simulations, or the barriers are so high that we have simply not encountered them yet for 
the size of triangulations which have been used in simulations of S 5 . 

In the following we assume that there are no problems with the Monte Carlo simulations 
for S* 4 related to the algorithmic unrecognizability. Since we are interested in studying 
the scaling relations near the critical point of the theory we are faced with another serious 

6 To be precise the phrase "combinatorial equivalence" means that there exists a triangulation which 
(up to graph isomorphism) is a common subdivision of the two triangulations. 

7 The way to search for such barriers is the following: Generate by Monte Carlo simulations a number of 
independent configurations for some large values of N4 . Now "shrink" (again by Monte Carlo simulations) 
these configurations to the minimal triangulation of S 4 , consisting of 5 4-simplexes. In case we never get 
seriously stuck in this shrinking procedure there can be no barrier separating two triangulations since we 
can first move to the minimal configuration and then out to another triangulation by the reverse set of 
moves. We have never observed that the triangulations get stuck in the process of a "reasonable shrinking 
procedure" . 
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practical problem when using only a standard set of the local moves. As will become clear 
from the results presented below, the typical geometry of the manifold depends strongly 
on the value of the gravitational coupling constant ki and becomes that of a branched 
polymer for large values of ki- The effect is a critical slowing down which makes the study 
of the branched polymer phase as well as the critical region very difficult. Local moves 
need many steps to change the geometric structure of the manifold and we observe very 
long autocorrelation times for quantities like (n(r)) N (33). To overcome this difficulty 
we enlarge the set of moves by global "baby universe surgery" moves which permit large 
changes of geometry and at the same time have large acceptance (in fact equal one). 
Similar moves were proposed earlier for two-dimensional systems [37] and led to a large 
reduction of the autocorrelation times. 

It should also be noted that the discussion of the algorithmic nonrecognizability deals 
only with the local moves. The theoretical implications of "baby universe surgery" are 
not known. 

4.1 Baby universe surgery 

The main new ingredient of our simulations is the "baby universe surgery" algorithm. 
The basic concept here is that of the "minimal bottleneck baby universe" introduced in 
Section 3. We recall that a "bottleneck" is a collection of five three-simplexes isomorphic to 
the boundary of a four-simplex, but not forming a boundary of any four-simplex belonging 
to the manifold. The "bottleneck" forms a boundary dividing the manifold in two parts. 
Since the original manifold has a topology of S* 4 both parts have a structure of S* 4 with 
the boundary ds^. The manifold can be split into two sub-manifolds with topology of S* 4 
by cutting it along the bottleneck and adding one four-simplex in each part to close the 
boundary ds^. This is the first step of the global "baby universe surgery" move. This 
step can be inverted: In each sub-manifold we remove a randomly chosen four-simplex, 
changing the topology of the sub-manifold to be that of S* 4 with a boundary ds^. The two 
parts can now be glued together and the boundary ds^ becomes a new "bottleneck" . The 
gluing can be performed in 5! ways. 

One can easily convince oneself that the process described above is self-dual: the 
numbers of simplexes and sub-simplexes remain unchanged. This means that there is no 
change in the action. The sizes of the "minbu" and the "parent" remain the same, only 
the position of the "bottleneck" is changed. The move also preserves the total number 
of "bottlenecks". Therefore the detailed balance condition is satisfied if the "bottleneck" 
from which we start the move is chosen randomly from the set of all "bottlenecks" of the 
manifold. 

We found it convenient to organize the "baby universe surgery" moves into a "global 
sweep". First we find all "bottlenecks" and form their list. The "baby universe surgery" 
is attempted for a randomly chosen "bottleneck" as many times as is the number of the 
"bottlenecks". The move is performed if the "minbu" size is bigger that a limiting value. 
In practice this limit was chosen to be 20. After the move the list of the "bottlenecks" is 
updated. The "global sweep" permits to obtain a distribution of the "minbus" for free, 
because this information is necessary anyway. 

The "baby universe surgery" moves are not ergodic. It is therefore necessary to sup- 
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plement them with the local moves. In the numerical study reported below the "global 
sweep" is followed by 9 standard sweeps. If the "global sweep" is not included, the au- 
tocorrelation times differ very much depending on the type of the measured quantity. 
Typically one has relatively short autocorrelation times for local quantities, like the aver- 
age curvature and very long autocorrelation times for the long-range quantities, like the 
average size of the system. Typically for a system with 16000 simplexes and ki ~ 1.2 (just 
below the pseudo-critical value) we get respectively 35 ± 5 and 2400 ± 200 sweeps. The 
autocorrelation times for the average size are in this case difficult to be measured, because 
a typical "time" dependence is so slow and in effect it can easily be underestimated in a 
short run. Autocorrelation times grow for larger systems as is usually the case (cf. the 
detailed discussion in the two-dimensional case presented in [37]) and for higher values of 
the coupling constant ki when the system enters the branched polymer phase. 

The "global sweeps" have small effect on local quantities, but dramatically reduce the 
autocorrelation time for the non-local quantities. For the system discussed above we get 
respectively 25 ±5 and 32 ±7 (the "global sweep" is treated as one of the sweeps here). The 
reduction of the autocorrelation time is even more spectacular in the branched polymer 
phase, where the measurements of the long-range quantities were practically possible only 
for small systems. The "baby universe surgery" algorithm permits us to study much larger 
systems (in our case up to 64000 four-simplexes) near the critical point. 

The important question when using the new algorithm is a comparison of the time used 
by the "global sweep" with that of the standard sweep. For a fixed value of the coupling 
constant ki the standard move, defined as N4 local moves, takes time proportional to 
N4. The performance of the "global sweep" depends on two parameters: the number of 
"minbus" and the average "minbu" size. Since the distribution of "minbus" is strongly 
peaked around small sizes the number of performed "baby universe surgeries" will depend 
on the cut-off in the accepted "minbu" size, which is a parameter of the program. In 
our calculations we found that the number of performed surgeries can easily be kept 
proportional to the system size N4. The average "minbu" size reflects the phase structure 
of the model and is a non-trivial function of N4 and ki (see fig. 2). It shows no N4 

1 /2 

dependence in the crumpled phase (small A^) an d seems to be proportional to N 4 in the 
so called branched polymer phase (in fact it provides additional evidence of the branched 
polymer interpretation as we will explain later). In effect the time used for a "global 
sweep" grows in this phase relatively faster than for a standard sweep. For ki = 1.6 this 
relative time increase is only 8/3 between 9K and 64K and even for the largest system the 
time of the "global sweep" remains of the same order as that of the standard sweep. 

4.2 Numerical setup 

As already mentioned above the details of Monte simulations of four-dimensional simplicial 
quantum gravity have already been reported elsewhere. The new ingredient in the present 
simulations, is the baby universe surgery algorithm combined with the increased theoretical 
understanding of scaling. The algorithm allows us to investigate in detail the elongated 
phase. We performed simulations for system sizes around 9K, 16K, 32K and 64K systems. 
By e.g. a 32k system we have in mind triangulations which are constrained to fluctuate in 
the neighborhood of 32.000 four-simplexes. We cannot keep the number of four-simplexes 
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Figure 2: The plot of the average minbu size versus ki for N4 = 9000 (triangles), 16000 
(squares), 32000 (pentagons) and 64000 (stars) simplexes. 



fixed in the simulations since the local moves in general change this number and they have 
to be included if the Monte Carlo algorithm is to be ergodic. We follow the method we 
used in our earlier work and enforce the fluctuations around the chosen size N% by adding 
to the action the extra term 

AS = A|A 4 - A 4 °| (48) 

with a small parameter A and choosing k 4 close to the critical value ^4(^2)- The value of 
A determines the width of the band around N® in which the number of four-simplexes is 
allowed to fluctuate. A priori it is not known if it is allowed to confine the variation of N4 
to a band close to the the desired average value for manifolds which are algorithmically 
recognizable, but the results reported in [28] support the idea that one can take quite a 
narrow band of a few thousand, or even less around the average value for all practical 
purposes. It is anyway possible to check that the results are numerically stable with 
respect to a given band-width. In practical calculations this width was taken to be 1000. 

In all the measurements we were monitoring the autocorrelation time and the typical 
run length was 200 to 300 r (r being the longest autocorrelation time). Near the pseudo- 
critical point it was sometimes necessary to increase this time to 500r. As was explained in 
the previous subsection the "baby universe surgery" algorithm permits us to measure the 
distribution of the "minbu" sizes at the same time when the "global sweep" is performed. 
Since the number of A4 fluctuates this distribution corresponds in fact to the distribution 
smeared over the band around N®. Since the band- width is small compared to A° we 
found that this smearing has no effect on the shape of the distribution and permitted us 
to get very good statistics for the baby universe distribution for free. 

The measurements of the distributions of (n(r)) N are the most time consuming part 
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of the program 8 . An accurate measurement requires averaging over the starting position, 
in principle over all four-simplexes of the manifold. In practice we found that we can get 
reasonable accuracy by restricting this averaging to a randomly chosen set containing 10% 
of the manifold four-simplexes. Below this fraction errors increase rapidly. Even with this 
reduction the measurement time for a manifold with size N4 grows like N%. For the 64K 
system 80% of the computer time is used by the measurement of {n(r)) N . 

4.3 The critical point 

The first question we want to address concerns the existence of a critical point k\ which 
separates two phases of simplicial quantum gravity in four dimensions. This transition has 
been known since the first simulations were performed. Recently it has been conjectured 
that it is a cross over and not a genuine transition [23]. Here we show that this is not the 
case. 

Define the average linear size of a universe of volume N4 in the following way: 



^) Ni = (^EE r(s 4 )n(r(s 4 )) ) , (49) 

where the average is performed over the different triangulations of S* 4 consisting of N4 
four-simplexes which are generated during the simulations. S4 denotes a four-simplex and 
n(r(s 4 )) denotes the number of four simplexes in the geodesic distance r(s 4 ) of S4. The first 
summation is over all four-simplexes, the next over all distances. For both summations 
we have a normalization factor I/N4. From the discussion in section 3 we expect in the 
case of a finite Hausdorff dimension that 9 

-c(rN 1/d *Y HndH ~ 1) 

WrJU-r^'e y 4 J , (50) 



JV 4 



i.e. we get: 



, d H /(d H -l) 



E^V- J ~KNl< d », (51) 



K= I dtt d *e- ctdHndH - 1 \ 



where the constant K is 



Eq. (51) is the "fixed 7V 4 " version of (N 4 ) r ~ r dn . 

In fig. 3 is shown a plot of the average radius of universes as a function of the coupling 
constant ki for various sizes of the universe (A^= 9K, 16K, 32K and 64K). It is clear that 
we see two phases, one where (r) N grows quite rapidly with N4 and one where there is 
almost no dependence on N4. It is also clear that we observe a drift towards higher values 
of &2 for the so called pseudo critical coupling /^(A^). One can define the pseudo critical 



8 They are performed only for N4 = N° , in contrast to the measurements of the minbu distributions 

9 Strictly speaking this formula is not valid: It is a hybrid of the short distance and long distance behavior 
which is convenient for deriving scaling relations like (51), but the power correction to the exponential 
decay is in general different from r dff_1 . 
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Figure 3: The plot of (r) Ni versus k 2 for 7V 4 = 9000 (triangles), 16000 (squares), 32000 
(pentagons) and 64000 (stars) simplexes. 



coupling in various ways. For spin systems one would use the peak in the specific heat or 
the peak in the spin susceptibility to identify the pseudo critical coupling. In such systems 
one expects a behavior: 

, x , f V~ x for first order transitions 
\P c {V)-p c \~\ , , 52 

| V fd tor higher order transitions 

where v is the critical exponent for the spin correlation function (not to be confused with 
v = 1/dh introduced above) and d the dimension of space. 

In the case of gravity it is natural to define the susceptibility as the second derivative 
of the partition function with respect to k 2 . In the same way as the magnetic susceptibility 
of a spin system can be given the interpretation as the integral of the spin-spin correlation 
function, we have seen that the second derivative of the free energy with respect to k 2 has 
an interpretation as the connected curvature-curvature correlation function: 

XR(k2, N 4 ) = dH ° gZ J^ NA) ~ (/ ^ 4 6vW6)vW6)) c , (53) 

where we have used a continuum notation. This susceptibility has indeed been used in 
the former studies. However, the peak is rather asymmetric for reasons to be explained 
below, and a precise location of the peak is not so clear. We find that the location of 
the pseudo-critical point k^N^) is more precisely determined from a measurement of the 
critical exponent j discussed in section 3. We will discuss these measurements in detail 
below. Here it is sufficient to show the best values of k 2 {N^) from these measurements in 
fig. 4. 
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Figure 4: The pseudo critical point ^2(^4)- 



We draw two conclusions from fig. 4: There is an infinite volume extrapolation k\ of 
^2(^4) an d it follows the finite size scaling formula (52): 

(54) 

The best fit is obtained for S = 0.47 ± 0.03 with the corresponding critical value k\ = 
1.336 ±0.006. It is tempting to conjecture that S = 1/2 as shown on fig. 4. The value of S 
indicates that the transition is not a first order transition. In case we use d = 4 since we 
consider four-dimensional manifolds 10 we get from (52) that the exponent v = 1/2. 

4.4 The entropy exponent 7 

We have already defined the entropy exponent 7 above as the subleading powerlike contri- 
bution to the partition function Z(k2]N^) or equivalently as the exponent which appears 
in the most singular part of the second derivative of the partition function with respect 
to the cosmological constant. A convenient way to measure 7 is by "minbu counting". 
Since we are already using a computer algorithm which identifies the baby universes this 
counting is for free, as already mentioned above. In the elongated phase ki > k\ there is 
a pronounced fractal structure and the distribution can be measured with great precision. 



10 It should be emphasized that precisely this aspect of finite size scaling is not yet clarified in the context 
of quantum gravity. For ordinary statistical systems the natural scale which enters in eq. (52) for higher 
order transitions is the linear dimension for the system. This linear dimension L is usually given to us by 
the geometry. Here we are integrating over geometries and only the total volume is unambiguously defined. 
This is why we reformulated (52) to be a function of the volume. The price is that only the combination 
vd appears. What value of d should we choose in the case of quantum gravity? The Hausdorff dimension 
du or the dimension d of the manifold. To our knowledge the answer is not known. 
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Figure 5: The measure values of 7 for &2 > ^2 f° r ^4 = 9000 (triangles), 16000 (squares), 
32000 (pentagons) and 64000 (stars). 



We find that 7 ~ 1/2. At this point it should be mentioned that 7 = 1/2 is a theoretical 
upper limit. The proof of this is identical to the corresponding one for two-dimensional 
quantum gravity coupled to matter [4]. In fig. 5 we have shown the values of 7 extracted 
from the measured distribution of baby universes. We simply fit the measured distribution 
to 



and make a lower cut such that only V4 > V^O) is included. Such a lower cut is expected 
since the theoretical formula is only valid for large V4. Typical cut-off value for V4 is 200. 
In fig. 5 we show the extracted value of 7 versus ki- It is seen that 7 ~ 1/2 for ki > k\ 
where k\ is approximately the value where the geometry changes according to fig. 3. In 
fact the measurement of 7 seems to give us the most precise signal for this change, a fact 
already utilized in fig. 4. Close to the transition 7 starts to drop and while the geometry in 
this region fluctuates quite a lot and makes precise measurements difficult, it seems that 
7 ~ at the transition point. Below the transition it is impossible to extract a 7 from 
the baby universe distribution 11 and as we will report below the asymptotic form (36) of 
the partition function is not valid in this region. It is possible to make a fit to (39) with 
the asymptotic form (37), and where a seems to change continuously from approximately 
1 to 1/4 when ki decreases from k\ to ki = 0. The precise determination of a was not 
possible, because acceptable fits are obtained for almost any value in this range. 

11 In an earlier letter [19] we reported about such measurements, but also remarked that the data for 
&2 < &f did not fit well to the straight line in the log-log plot corresponding to eq. (55) for large V4. 
We attributed it to the bad statistics of large baby universes, which are sparse in this phase. However, 
the evidence now points to the non-existence of a 7 in this phase and in fact a distribution like (36) is 
disfavored for &2 < k%. We therefore retract the claims in [19] for &2 < k% concerning 7. 



log Af(V 4 ) = (7 - 2) log[V 4 (l - V 4 /N 4 )] + const. 



(55) 
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Figure 6: (r) N /y/Wl versus k 2 for N4 = 9000 (triangles, 16000 (squares), 32000 (pen- 
tagons) and 64000 (stars). For k 2 = 1.4 and 1.6, i.e. away from the critical point 
k 2 = 1.336, the values are independent of N4. 



In the next subsection we will provide additional evidence that 7=1/2 for k 2 > k 2 . 
4.5 The distribution (n(r)) 

The measurement of (n(r)) N is easy in the numerical simulations. One simply counts the 
number of four-simplexes in shells of distance r, and one can recursively work one's way 
out from r to r + 1 until the whole manifold is covered. Taking the average over different 
starting positions and the average over independent triangulations with a fixed N4 will 
provide us with (n(r)) N . 

As discussed in section 3 we expect from theoretical considerations that if the ensemble 
of piecewise linear manifolds has a characteristic internal Hausdorff dimension the function 
{n(r)) N should obey the scaling law: 

n(r,N 4 ) = ±-( n (r)) Ni = Np r^r/A^). (56) 

In addition we expect that the short distance behavior of n(x) to be like x dH ~ 1 and the 
long distance behavior to be like exp( — c x dH ^ dH ~ 1 ^) . 

In fig. 6 we have tested the scaling behavior of n(x) for {r) N and k 2 > k\. If the form 

(56) is satisfied (r) N = J2 r rn ( r i N4) ~ N^ ds . In the figure djj = 2 and it is clear that 
we have perfect scaling for k 2 not too close to /^(A^), i- e - for k 2 above k 2 = 1.34 there 
should be scaling for all values of N4. 

The data are sufficiently good that we can in fact fit to a detailed functional form for 
n(x). For k 2 = 1.29 and N4 = 64.000 we show in fig. 7 the data and the function: 

n(x) ~ xe~ cx2 . (57) 
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The data points are the dots (the error bars are less than the dots) and the curve is just 
n(x) from (57) with only one adjustable parameter (c). It is seen that the agreement is 
almost perfect in spite of the fact that ki = 1.29 is just above the transition point (the 
pseudo critical point for N4 = 64000). If we move further into the elongated phase the 
agreement gets even better. This implies in particular that the distributions for the same 
&2 can be scaled on top of each other by using x = r/yW^ in this phase. In fact the 
agreement is so good that it would be impossible to distinguish systems of different sizes 
for &2 = 1.4 and 1.6. 

The conclusion so far is that the elongated phase ki > k\ is characterized by 7 = 1/2 
and dfj = 2 and that 

G(r; N A ) ~ N^' 2 re^ 2 ^ (58) 

From this we can construct the two-point function: 

00 

G(r;k 4 ) = J2 G(r;N 4 )e- k * N * ~ e - c V^>\ (59) 

JV 4 = 1 

This two-point function is precisely the two-point function of a theory of branched polymers 
[4] and the picture which emerges is that the typical manifolds in this region of the coupling 
constant space are like small tubes which can join and branch in tree-like structures. 

The transition away from the branched polymer region takes place increasingly abruptly 
with increasing N4 and we conjecture that the phase in the infinite volume limit continues 
all the way down to k\. In fact this picture would be consistent with the observed finite 
size scaling exponent v = 1/2 introduced below eq. (54), if we identify this exponent with 
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the mass exponent of the two-point function which is v = 1/djj- At the transition point k 2 
there could be a different critical behavior. It is quite difficult to perform measurements 
exactly at the transition. The geometry fluctuates a lot, but as already mentioned the 
minbu counting suggests 7 ~ 0. Measurements of dfj suggests dfj ~ 4, but we must insist 
that the data do not allow a clear extraction of dfj. 

The measurements of (n(r)) N for k 2 < k 2 show completely different characteristics than 
for k 2 > k\. For large r there is no N4 dependence at all and the function seems to fall 
off exponentially: 

(n(r)) Nt ~e- m W r . (60) 

As discussed above such a behavior is consistent with the notation of infinite Hausdorff 
dimension. However, if this notation should be "genuine" it should also reflect itself at 
small distances, i.e. contrary to the situation for smooth four-dimensional manifolds with 
negative curvature it should not be possible to have {n(r)) N = r 3 for small r. We have not 
been able to fit to any sensible power dependence r dH ~ 1 for small r. A plot of log (n(r)) N 
shows indeed a linearly growing function of r up to some ^(A^) which can be identified 
with both the peak of the distribution (n(r)) N and with the average value (r) N . ro(N4) 
or {r) N shows only a very weak dependence of A4 and a reasonable fit to {r) N is 

(r) Ni = a{k 2 ) + b{k 2 ) log N 4 . (61) 

This again gives some support to the idea that the Hausdorff dimension is infinite in 
this phase. Recall from section 3 that an infinite Hausdorff dimension will result in an 
exponential dependence like in (60) if it appears as a limit of distributions with finite 
Hausdorff dimension. Similarly eq. (61) would appear as the limit of (^)jy ~ 

]yi/dH f or 

dfj —7- 00. 

The observation that (n(r)) N grows exponentially from r ~ 6 out to r ~ ro and then 
falls off exponentially indicates that we deal with an infinite Hausdorff dimension at all 
distances and it is easy to get a quite good "phenomenological" fit to (n(r)) N which 
incorporates both these features by choosing e.g.: 

(n{r)) Ni ~ exp (-m x {k 2 )r - c 2 e~ m ^ r ) . (62) 
It will grow like e^ C2 ' m2 ~' mi ^ r for small distances and fall off like 

M r ))iV 4 ~ c i e~ mi ^ r - C2 e-( mi ( fc2 ) +m2 ( fc2 ) r + • • • (63) 

for large distances, while a A4 dependence in the coefficient c 2 would explain the observed 
A4 dependence of ro- 

The data and a fit of the form (62) are shown in fig. 8 for k 2 = 1.26 and A4 = 64000, 
i.e. right below the transition to the crumpled phase, where the fit is worst. But even so 
close to critical point (62) works quite well over the whole range of r > 6. It should be 
mentioned that the coefficient in front of the second exponential in eq. (63) is negative. 
This implies that the term cannot be given the interpretation as an additional heavier mass 
excitation. However, just looking at the long distance tail the distribution (r) N allow us 
to determine nil, m 2 and c 2 from (63). On the other hand we can determine c 2 m 2 — mi 



23 




0.05F 



20 40 60 



80 



100 



120 



140 



r 



Figure 8: exponential fit (curve) of the form (62) to the measured n(r; N4) (dots, error 
bars less than dot-size) in the crumpled phase (N4 = 64000, k 2 = 1.26). 



from the short distance exponential growth alone and find good agreement. This indicates 
that long and short distance behavior are intertwined in the case of infinite Hausdorff 
dimension, as they are in the case of finite Hausdorff dimension, where dfj appears both 
in the short distance and long distance expression for (r) N (see e.g. (50)). 

We conjecture that the internal Hausdorff dimension is infinite for k 2 < k 2 .. 

It is worth to emphasize that (60) is inconsistent with even a logarithmically weak 
scaling of the mass of the two-point function G(r; k^) for A&4 — > 0. As mentioned in 
section 3 it would require at least a logarithmic dependence like e _r / log ( JV / r ). Our data 
seem not to favor such dependence and this implies that the two-point function (by the 
usual discrete Laplace transformation) will be: 



The mi(k 2 y& do not scale to zero as k^ — > Kl(k 2 ), but in fig. 9 we have shown their limiting 
value at k1(k 2 ) as a function of k 2 — /^(A^). The data seem to fall on universal curves 
and the lowest mass m\{k 2 ) scales to zero as k 2 — > k\. If we extract a critical exponent 
v (not to be confused with the former introduced mass exponents v = 1/dh and v from 
finite size scaling) for mxik^) we get: 



The fact that we have a mass excitation which scales to zero at the critical point k\ is 
a strong indication that the system will develop genuine extension and a finite Hausdorff 
dimension at the critical point. 



G(r; k 2l k - 4) ~ A x e 



mi (k 2 ) r 




(64) 




(65) 
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Figure 9: The behavior of the masses mi^) and mi^) + ^2(^2) from (62) in the 
crumpled phase as a function of k\ — ki- 



4.6 Miscellaneous 

Let us collect here some of the observations which were scattered in the presentation of 

the data, but which did not fit into the discussion earlier. 

The first observation is that in the elongated phase the average "minbu" size grows as 
1/2 

N 4 , N4 being the total volume of the universe. This was mentioned above and is shown 
in fig. 10. Is it understandable? The answer is yes if we assume that 7 = 1/2. If we assume 
there exists a 7 at all we calculated earlier the probability distribution for minbu universes 
of size V4 to be N4 [V^(l — V4/./V4)] 7-2 (see eq. (40)). The average size is therefore 

JV4/2 1/2 
<^4>jv 4 = £ y 4 7V 4 [V 4 (l - V A /N A )\<- 2 « N2 / dW~\l - ty~ 2 (66) 
V i = l Jo 

if we assume 7 > 0. The data support the value 7 = 1/2. 

The next observation concerns the behavior of the pseudo-critical point /^(A^) as a 
function of N4. According to eq. (54) it approaches the true critical point k\ as l/y/N^ 
and we conjectured that it corresponded to critical exponent v = 1/2. Do we have such 
mass exponents around as a function of ki ? The answer is yes: The mass raxik^) in 
the crumpled phase seems go to zero as \ki — k^l" , where the exponent v ~ 1/2. A 
closer analysis of the data of for instance the average curvature and the susceptibility 
X-r(^2j N4) indicates that the major part of the finite size effect comes from the crumpled 
phase. The curve corresponding to X-r(^2j N4) shows a clear asymmetry as a function of 
ki around the peak at /^(A^). It grows relatively slowly from the crumpled phase when 
^2 < ^2(^4) approaches /^(A^) an( l drops abruptly for ki > /^(A^). The asymmetry 
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Figure 10: The scaled average minbu size (b) N denned in eq. (66) versus ki- 



can now be partly understood: In the crumpled phase we have no masses scaling to zero 
as &4 — > ^4(^2)- Only when ki approaches k\ will the lowest mass scale to zero. In 
the elongated phase we always have a "mass" excitation to (£14) which scales to zero as 
&4 —7- ^4(^2) like (&4 — k\(k2)) l l 2 ■ This true for all values of ki > k\. It is natural to 
conjecture that the finite size dependence of the pseudo-critical point ^2(^4) is monitored 
by the mass raxik^) in the crumpled phase. 

5 Discussion 

5.1 Summary of observations 

Four dimensional simplicial quantum gravity as defined in this article depends on two bare 
coupling constants, ki and k^. The first serves as the inverse of the bare gravitational cou- 
pling constant while the second one in the formal continuum limit is a linear combination 
of the inverse gravitational coupling constant and the cosmological constant. We expect at 
least &4 to have an additive renormalization, i.e. the continuum limit of the theory should 
be defined as a function of k^ — k\{k2). In the (^21^4) coupling constant plane there is a 
critical line k^ = ^4(^2) where such a limit can be taken. This leads to an infinite (lattice) 
volume of the theory. The renormalized cosmological constant will be defined as a function 
of &4 — k\{k2). It might not be possible to define an interesting continuum field theory 
for all values of ki- Recall the analogous situation for a spin model, e.g. the Ising model, 
where there is only a single point where a continuum limit can be taken. 

We have observed two distinct phases of the discretized theory. One phase, denoted 
the elongated phase, was characterized by specific scaling properties of the two-point 
function, which allowed us to conjecture that the Hausdorff dimension of the ensemble 
of (piecewise linear) four-manifolds considered was dfj = 2 while the entropy exponent 
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7, which determines the proliferation into baby universes, was 7 = 1/2. The detailed 
investigation of the geometry in this phase showed perfect agreement with a model of 
fractal structures denoted branched polymers. This phase was present as long as the 
inverse bare gravitational coupling constant ki > k^, i.e. for small ("semiclassical") values 
of the bare gravitational coupling constant. 

For &2 < k\ we nave a completely different phase. From the numerical data we were 
led to conjecture that the Hausdorff dimension dfj is infinite and that there is no entropy 
exponent 7, in agreement with recent results obtained for A; 2 = [28,29]. In this phase we 
can say that there is effectively no metric associated with the volume: The volume enclosed 
in a geodesic ball of radius r grows exponentially with r, but contrary to the situation 
for smooth spaces with negative curvature this is true down to the smallest distances of 
r ~ 3 — 4. This implies that the main part of a typical universe can be reached in few 
steps, and we can view the universe to be of Planck size. This is possible because there 
exist vertices of very high order in the typical triangulation, and such vertices ensure that 
a large part of the universe can be connected in few steps. Nevertheless the exponential 
decay of the two-point function allows us to identify a mass which scales to zero when 
^2 — > &2- This strongly suggests that the ensemble of four-manifolds will develop a genuine 
extension when ki — > k\. 

It would be very interesting if we could determine the precise value of the Hausdorff 
dimension at the critical point k^, as well as the entropy exponent 7 if it exists in this 
point. The present data hint that 7(^2) ~ but it is clear from the curves of fig. reffig4 
that the correct should rather be: in the critical region, where 7 changes very fast and 
probably even jump from 7 = 1/2, we have measured values of 7 ~ 0, but never values 
below zero. This does not rule out that increased fine tuning of ki — > k\ from above would 
allow us to reach negative values of 7. 

5.2 Comparison with Regge discretization 

Apart from the dynamical triangulation approach to quantum gravity, there exists another 
approach which here we shall denote Regge quantum gravity (RQG). The approach is 
completely different in philosophy from simplicial gravity. While simplicial gravity should 
be viewed as a discretization directly in the space of equivalence classes of metrics, RQG 
works with a fixed triangulation and varies the length of the links. In this way many 
configurations might correspond to the same metric and in principle a Jacobian should 
be included in the functional integral. This may be a serious problem for the approach. 
On the other hand it has the advantage that it allows for a weak coupling expansion 
which can be important if we want to relate to a continuum formulation. For an review 
we refer to [41]. Recently some problems have been reported both in two-dimensional 
gravity where the simplest choices of measures do not reproduce the known continuum 
results [38,39] and in four dimensions, where a drastic dependence on the choice of the 
fixed triangulation used in the computer simulations was observed [40]. If we ignore these 
problems and just refer to the results as reported in [41] there seems to be some similarity 
in the results obtained from RQG and the ones reported here. First of all there are two 
phases in both approaches. The elongated phase is not particularly well defined in RQG, 
but it still has some flavor of the elongated phase seen in simplicial quantum gravity. 
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In the other phase a mass gap is observed, precisely as in our simulations and the mass 
scales to zero with approximately the same exponent as v in eq. (65). There seem to be 
divergent opinions concerning the sign of the volume- volume correlation function. In [42] 
it is reported negative, while in [43] positive. In simplicial gravity the correlation function 
is by definition positive. 

There are also differences in the two approaches. In simplicial quantum gravity the 
naively defined average scalar curvature does not scale to zero at the critical point. It is 
not presently known if this is due to the fact that the curvature operator is a composite 
operator which requires a subtraction before contact can be made with continuum. In 
that case the average value of the curvature should be absorbed in a redefinition of the 
cosmological constant. Whatever is the explanation, it disagrees with the measurements 
from RQG where the average curvature is seen to scale to zero when ki — > k\. Another, 
more drastic difference is indicated by the Hausdorff dimension in the strong coupling 
phase (i.e. the crumpled phase in simplicial quantum gravity). As reported above it 
seems to be infinite in crumpled phase. In RQG it is reported to be close to four. 

5.3 Continuum interpretation 

As we have shown above simplicial quantum gravity is a perfectly well defined statistical 
theory where it is possible to define and discuss the continuum limit and derive scaling re- 
lations in a way similar to what is done in the ordinary theory of critical phenomena. Here, 
as in statistical mechanics, it might be a problem to identify a continuum theory which 
corresponds to the observed critical phenomena. The problem is that we are presumably 
far away from any "naive" continuum limit. This is not necessary an undesirable situa- 
tion. Any "naive" continuum limit of Euclidean quantum gravity is a disaster because the 
action is unbounded from below. A number of cures have been suggested: The addition 
of suitable higher derivative terms can bound the action and at the same time make the 
theory renormalizable [44]. A sophisticated choice of functional measure starting from the 
formal expression in spaces with Minkowskian signature can also tame the unbounded- 
ness [45] (see [46] for a recent review). Finally the problems with the unboundedness of 
the action should disappear in a topological phase of quantum gravity, where the metric 
should play no role. Three-dimensional quantum gravity in the formulation of Witten [47] 
provides an instructive example of how a "minor" change from metric to drei-beins allows 
us to view the theory as topological and how the rotation to Euclidean signature results 
in a bounded action. Further the Euclidean theory has a beautiful discretized realization, 
the so called Turaev-Viro theory [48], which in spirit is not far from simplicial quantum 
gravity. In fact it is presently not ruled out that three-dimensional quantum gravity in 
the crumpled phase can be given a representation as topological gravity. 

In simplicial quantum gravity we have no problems with an unbounded action as long 
as the lattice volume is finite. However, part of the regularization is performed directly 
in the equivalence classes of metrics. This makes it difficult to decide which of the above 
scenarios is realized by the regularization, and it could depend on the choice of the bare 
(inverse) gravitational coupling constant ki- Due to the lack of analytic solutions of the 
discretized theory as in two dimensions, we have to take the numerical data and look for 
scaling behavior which can be associated with one of the above continuum models. 
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In the elongated phase we understand now in detail the statistical theory, namely the 
one of branched polymers. It can be viewed as the ultimate dominance of volume-volume 
fluctuations and therefore of the conformal factor. This makes it natural to compare with 
the models of four dimensional quantum gravity which study the conformal sector based 
on the trace anomaly [49,50]. This approach, which we denote conformal four-dimensional 
quantum gravity (CFQG), is inspired by the corresponding treatment in two dimensions, 
but the same is true for our generalization of the successful two-dimensional simplicial 
quantum gravity theory. CFQG studies the fluctuations of the conformal factor 
space-times of the form: 

9nv = e 2 °( x %, v {x), (67) 



where g llv {x) is a fiducial metric. The theory has an infrared fixed point and at this fixed 
point the field a receives quantum dressing. The effective action can be written as [51] 
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In this formula the first term is the trace anomaly induced action and k = 87r£r, while the 
second term is the Einstein action action, a denotes the quantum gravitational dressing 
coefficient, A4 a fourth order differential operator with respect to the fiducial metric g, 
G the Gauss-Bonnet term, while Q 2 is a kind of analogue of the two-dimensional central 
charge. Finally f(x) = x(l + Ax + 6a; 2 ). In addition the partition function at fixed volume 
will scale like [51]: 

z(k,v) ~ v^zLv 1 ' 2 } 
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The dressing coefficient a is related to the gravitational "central charge" by: 



a 
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In particular we expect from (69) that the approach to the infinite volume limit is: 



k ~ 



yl/2' 



i.e. 



1/2' 



(71) 



where the last equation is the translation to the discretized notation of simplicial quantum 
gravity and where we have identified the critical point k\ observed in simplicial quantum 
gravity with the infrared fixed point of CFQG. 

We are now in a position to compare with the "numerical observations" reported 
above. At the critical point we found l&^A^i) — k?,\ ~ , where 5 = 0.47 ± 0.03. 
The finite size dependence of the pseudo critical point is derived from finite size scaling, 
where the assumption is that the (singular part with respect to A^) of the free energy 
F(k2, N4) = log Z(k2, N4) can be written as a function F(N^\k2 — k^l). A glance at (71) or 
(69) shows that the measured S indeed agrees with the prediction of CFQG. Unfortunately 
it is difficult to measure 7(^2) ver y precisely at the critical point. Our measurements 
indicated 7^0 which corresponds to a value Q 2 = 16/3 and a = 4/3. According to 
the hypothesis of "infrared conformal dominance" [51] this value of Q 2 should arise by 
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integrating out the transverse gravity mode, thereby creating an effective action of the 
same form as (68) when expanded in powers of derivatives. 

However, there is some ambiguity in this comparison since it is possible to add an 
independent volume term to (68) [46,52]: 

S eII ^S eII + A J rf 4 vf e 4 ^. (72) 

In this case the scaling is changed to 

Z( K ,V)~V^ 3 z( K V a W), 7 = 2-^, (73) 

where 

o = 1 - >/l - 8/Q 2 4/Q 2 

' 4/Q 2 ' a 2/Q 2 

and we can eliminate the unknown Q 2 : 

' = §-ip-7-V> + (>-T)»). P5) 

If we compare (73) with the observed finite size scaling it follows that we only obtain the 
observed a/2/3 = 1/2 if j — > — oo at k\. As remarked above this can not be ruled out 
from our measurements. This would correspond to the "semiclassical" region of CFQG 
where Q 2 — > oo and the analogous regime in two-dimensional quantum gravity would be 
c —7- — oo. If we take 7 ~ we get a value a/2/3 = 0.29 which is below the measured value. 
It is very unfortunate that we presently are unable to determine 7 more precisely at the 
critical point. 

When we move into the elongated phase of simplicial quantum gravity it seems that 
we have reached an extreme limit. In the framework of extended CFQG defined by (72) 
one would be tempted to view this as the transition where 7—7-0 and /3 ceases to be real 
for Q 2 —7- 8. This transition has a flavor of the c = 1 barrier in two-dimensional quantum 
gravity and it has often be conjectured that c = 1 (where 7 = in two-dimensional gravity) 
marks a transition to branched polymers. 

To summarize: The comparison of simplicial quantum gravity and CFQG is not en- 
tirely unambiguous, partly due to the numerical uncertainty associated with the location 
of the critical point and partly due to the approximation which is inherent in CFQG by 
assuming that all effects of the transverse degrees of freedom can be absorbed into Q. 
However the qualitative agreement is very encouraging. 

In the crumpled phase the situation is more controversial. Seemingly the Hausdorff 
dimension is infinite (or at least quite large). In addition it has not been possible to define 
an entropy exponent 7 and formally 7 = — 00. This is in agreement with "observations". 
In the crumpled phase we observe only few and small baby universes. If we move deep 
into the crumpled phase we find in addition that the naively defined curvature is negative. 
Does this imply that we observe smooth manifolds with negative curvature? On such 
manifolds one would expect the volume to grow exponentially with the radius of geodesic 
balls. The answer seems to be no. Rather the scenario is one where a few vertices have 
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very high order and are linked to almost all other vertices in the manifold. On such 
piecewise linear manifold the link distance between any two vertices will be only a few 
lattice units 12 and we see only a very weak dependence of the average radius on the total 
volume. Nevertheless we know from theoretical considerations that we can define a mass 
by the exponential decay of the two-point function and from the numerical simulations 
there are indications of a hierarchy of masses. It is tempting to conjecture that the concept 
of volume and consequently of metrics cannot play a major role in this phase, even in spite 
of the ability to define a two-point function. In this case we approach a situation where 
the linear extension of any universe of finite volume is of Planckian size and where a more 
adequate description might lead to a theory which has reference only to the topology 
of the underlying manifold. It is interesting to note that if we choose ki = 0, i.e. we 
perform a pure counting of manifolds with no additional weight, we are deeply into the 
crumpled phase: the generic piecewise linear manifold is of Planck scale. As we move 
towards the semiclassical region (the bare gravitational coupling constant 1/&2 ~~ ^ 0) the 
ensemble of manifolds undergoes a phase transition: Above the transition the ensemble 
has a finite Hausdorff dimension dfj. In fact dfj = 2. It is of course a most interesting 
question to determine the Hausdorff dimension at the transition point. This point seems 
a candidate for defining a continuum non-perturbative theory of quantum gravity with 
a finite Hausdorff dimension larger than two and it has the tantalizing interpretation as 
transition point between the entropically favored Planckian size quantum universe and the 
two-dimensional fractal universe dominated by the conformal mode. 

5.4 Outlook 

It is worth to recall that the action used in simplicial quantum gravity is extremely simple. 
We consider this as a virtue of the theory. It has certainly been very useful in the studies 
of two-dimensional quantum gravity, and there is hope that it will allow us to solve the 
four-dimensional theory analytically. Until then numerical simulations can be quite useful, 
serving as inspiration for analytic approaches, and powerful enough to give us important 
hints about the solution of the theory. The indications that the elongated phase is a 
branched polymer phase 13 and that the crumpled phase might correspond to a topological 
phase corroborate on the hope that it might be possible to solve the model analytically. 

As already emphasized above it seems very important to develop better algorithms 
in order to study in detail the neighbourhood of transition point. The baby universe 
surgery algorithm has been crucial for this numerical study, but there is still room for 
improvement. 

A most interesting question which we have not touched upon is that of topology 
changes. Hopefully the full power of dynamical triangulations will show us a way to 
deal with these changes in a quantitative way, as happened in two-dimensional quantum 
gravity. 

12 The corresponding "four-simplex distance" is in general somewhat larger, but it exhibits the same 
scaling dependence as the link distance. 

13 It should be emphasized that the branched polymer phase is much clearer in four-dimensional space 
than in two-dimensional gravity coupled to matter. 
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